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We investigate simplified models of dark matter with scalar mediators at hadron colliders using the 
final state topology with 2 jets and missing energy. These models can arise in a wide variety of 
BSM scenarios including the possibility of the mediator mixing with the Higgs. Our aim is first to 
determine the projected reach of the LHC and the future circular hadron collider for excluding such 
models, and we also compare these to the relic density and direct detection constraints. We use the 
kinematic distributions to extract information on mediator masses at colliders. At the 13 TeV LHC 
we can probe mediator masses up to 750 GeV, and at a 100 TeV collider the reach is increased to 
2.5 TeV mediators. We also explain how individual models with different values of mediator masses 
can be differentiated from each other. 


I. INTRODUCTION 

There is overwhelming cosmological evidence for the existence of dark matter (DM) with a density about five 
times larger than that of ordinary matter. But despite a lot of ongoing effort in direct and indirect detection, and 
collider experiments we still do not know the fundamental nature and composition of DM. None of the particles in 
the Standard Model can provide a dark matter candidate while many beyond the Standard Model (BSM) scenarios 
do. A common feature of many BSM extensions are the prediction of Weakly Interacting Massive Particles (WIMPs): 
if a particle has weak scale mass and interactions, the abundance after standard thermal freeze-out would be close to 
the observed abundance of dark matter. Specific realisations of DM candidates in different BSM models can have a 
rich and interesting phenomenology. However, assuming dark matter can be produced at colliders, a generic feature 
is the existence of signatures with missing energy, as stable dark matter particles leave the detector unobserved. 

To discover missing energy signals at colliders, one recoils the invisibly decaying particle against reconstructable 
objects. In the simplest and most direct case, this is a visible mono-object, and such searches for mono-jets and 
mono-photons have been carried out at Run 1 of the LHC m- These studies have so far not discovered any evidence 
for an excess of missing energy events and can in parts of the parameter space be as or more constraining than limits 
from direct and indirect detection HD- It is thus important to formulate and extend the searches of dark matter at 
the LHC in Run 2 and beyond. 

Dark matter can be produced at colliders via an exchange of a mediator particle which connects the colliding 
SM partons to the dark sector. A viable and simple approach to characterise and interpret dark matter searches 
at colliders relies on using simplified models with four basic types of mediators: vectors, axial-vectors, scalars and 
pseudo-scalars (see white papers [8l|9] for early reviews and references). The mediator is a dynamical degree of freedom 
in this approach, and this is the correct description for dark matter searches at the LHC as the energy transfer in 
the collision can typically exceed mediator masses. Following the Higgs discovery, there is a renewed interest in the 
role of scalar degrees of freedom and the possibilities provided by extended Higgs sectors in searches for new physics. 
Of particular interest to dark matter searches are the models with scalar and pseudo-scalar mediators whose reach 
at the LHC was studied recently in [inHH]. It was found that the LHC at 14 TeV will provide a complementary 
coverage to the low-energy experiments in dark matter searches, and can be the only experiment to probe dark sectors 
if the invisible particles produced are not stable at cosmological scales. These studies have been performed using the 
so-called mono-jet topology m- 

In this paper we will study simplified models with scalar mediators in the 2-jets plus missing energy topology to 
determine their collider limits and the discovery potential by analysing the kinematics of the final state jets. For 
scalar mediators mono-jet searches are predominantly relying on gluon fusion production mm- Now the presence 
of a second jet allows for a more non-trivial kinematic distribution to characterise the final states. Hence the VBF 


* Electronic address: valya.khoze@durham.ac.uk 

^Electronic address: g.o.i.ro@dur.ac.uk 

^Electronic address: michael.spannowsky@durham.ac.uk 



2 


cuts can be imposed which suppress the gluon fusion production channel for scalars. This makes the weak boson 
fusion processes dominant instead and allows to capture mediators with suppressed couplings to fermions. Thus, the 
kinematic information in the 2 jets+MET final state should allow to probe better the mediator masses and also to 
give a handle on their interactions with electroweak gauge bosons. In a slightly different context, the idea of exploring 
the two jets kinematics to learn more about the SM~DM interactions has also been implemented in m- 

We consider a simplified model with a scalar mediator whose SM couplings proportional to the SM Higgs. We start 
by defining the models in Section El before commenting on non-collider dark matter phenomenology in Section m 
Then, we move on to consider the phenomenology at the LHC in Section m and at a future 100 TeV collider in 
Section 


II. MODELS 


In unitary gauge the Standard Model (SM) contains just a single scalar-field degree of freedom, the neutral scalar 
Higgs h. At tree level h interacts with the massive vector bosons, and and all the SM fermions. The linear 
interactions in h of the Higgs boson with other SM particles can be written in the form. 


rSM ^ 


w+w->^ + ^ — ff 

V ^ V V 

f 


h . 


( 1 ) 


We want to extend the SM by introducing a scalar mediator particle (p which couples to the SM degrees of freedom 
as well as to fermionic dark matter particles x via. 


A -5DMXX<(>- 


( 2 ) 


For the purpose of this article the spin of the dark matter particle is not relevant, i.e. the dark matter particle 
could be instead a vector or scalar particle. There are two types of settings where the additional scalar p can appear 
in interactions with the Standard Model. First, it can be an additional Higgs doublet, for example coming from a 
two Higgs doublet model, or more generally any scalar field transforming non-trivially under the SU{2)j^ of the SM. 
Alternatively, the (p scalar mediator can be a singlet under the Standard Model. In the latter case, it interacts with 
the SM degrees of freedom only via the mixing with the SM Higgs h. In this case, the interactions of p with the 
SM are subject to experimental constraints on the mixing angle sin^ 9 < 0.15 (see mun]) arising from experimental 
bounds on the SM Higgs to invisibles decays and other Higgs data. 

First in section IIA we consider the more constrained singlet-mixing case, and next in section IIB| we define the 
less-constrained generic Higgs-like scenario. The upshot is that both of these cases will be described by the same 
simplified model of Eq. (12) with the scaling parameter k being either unconstrained k ^ 1 or small n < 0.15. 


A. The singlet mixing model 

Here p is a. Standard-Model singlet neutral scalar. The visible SM sector and the ‘invisible’ y sector are coupled to 
each other only via the mixing between the two neutral scalars, p and h, as in the Higgs portal model. The states of 
definite masses, hi and / 12 , are 


h = hi cos 9 + h2Sm9 , p = —hi sin 9 -|- /i 2 C 0 s 6 *, (3) 

where 9 is the mixing angle. Combining Eqs. 0-@ we obtain a Simplified Model for invisible Higgs decays involving 
two Higgs-like neutral scalars hi and / 12 : 




hi,h2 
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(4) 


The first scalar mass eigenstate hi plays the role of the observed SM Higgs boson and we also assume that the 
mediator /12 is always heavier than the SM Higgs, 


> mhi = 125 GeV. 


( 5 ) 
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With this Lagrangian we can produce ^2 as in the SM via both gluon fusion and vector boson fusion mechanisms, 
with the corresponding SM cross-sections rescaled by sin^ 9. Similarly the 125 GeV Higgs scalar hi production rates 
are rescaled relative to the SM by the factor of cos^ 9 which is ~ 1 for sufficiently small values of the mixing angle. 

If both mediators can be produced on-shell, in either channel the cross-section for xx-l-2 jet production in the 
narrow width approximation can be written as. 


where ah^ is the production cross-section for hi -\- 2 jets and Jivh-^xx branching ratios. 


where 
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(9) 


For dark matter mass below the kinematic threshold of both mediators, 2my^ < mh^ < mh 2 both mediators can be 
on-shell and, in principle, both channels for the dark matter production are open, but the lighter Higgs will dominate, 
as can be seen from ([^-(|^, 
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( 10 ) 


The SM Higgs has a very narrow width of 0.0068 GeV and due to the limits on the Higgs to invisible branching ratio 
we know that this width, can not increase by more than 35%. The reason the Higgs width is so small is due to 
the fact that all the fermions are coupled to the Higgs via Yukawa couplings so we cannot have light fermion with 
large coupling. The total decay width of the second scalar, on the other hand can easily be large as can be 

inferred from Fig. Even for ^dm =0.1 the total width of /12 will be an order of magnitude larger than the hi Higgs 
width. Hence, for light dark matter, when both channels are open, only the hi Higgs mediator is relevant. 



FIG. 1: The decay width of h 2 into XX with ^dm = 1. 


For heavier dark matter, mh^ < 2.m-^ < only fh® ^2 channel is open and it is efficiently described by the 

simplihed model. 
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Finally, if dark matter masses are higher than it cannot be produced via an on-shell mediator exchange, and the 
resulting rate of its production is too small to be observed. 

Current limits on sin^ 9 mainly come from two sources, the Higgs signal strengths and the electroweak precision 
tests, for recent papers see [HHiH]. Limits from Higgs signal strength measurements constrain cos^ 9 directly m- 
This leads to a bound sm9 < 0.44, independent of the mass of li 2 - The electroweak precision tests, mainly the W 
boson mass give a mass-dependent constraint on sin 9 shown in figure 3 in Ref. [15j . In the mass range around 1 TeV 
the limit becomes sin 9 < 0.3. We also note that the limits coming from a non-observation of the second SM-Higgs-like 
state are not directly applicable for /i 2 in our case due to its large branching ratio to invisibles. 

We will only consider these limits in the context of the singlet-mixing simplified model with the K-parameter 
K = sin^ 9. In the simplified model framework we do not know what other particle content there is and additional 
degrees of freedom could still modify both the SM Higgs signal and the loop corrections to the W-mass. 

Recent discussion of theory models for dark matter based on mass mixing between the scalar singlet mediator and 
the SM Higgs can be found in Refs. [I9U24] . 


B. Generic Higgs-like scalar mediator model 


More generally, scalar mediators to dark sector can also arise from an independent additional Higgs doublet or 
Higgs multiplet, for example in the two-Higgs-doublet models. We choose the simplified model for a generic scalar 
mediator by assuming that it has the same interactions with the SM vector bosons and fermions as the SM Higgs, 
but scaled by an overall scaling factor k which is a free parameter of the simplified model. 


£ = I + 


Ml 




- ffDMXX'/' - - m^xX- 


( 12 ) 


In general, the scalar mediator can couple with a different strength to the SM vector bosons and to SM fermions, thus 
introducing additional parameters into the simplified model (12|. For clarity and simplicity we will use the minimal 
model (12) with a single scaling factor. Here k = 1 corresponds to the normal SM Higgs couplings. In general 


we consider values of k ^ 1 since it is difficult from a model-building perspective to increase the coupling to gauge 
bosons with additional Higgs singlets or doublets. The simplified model for the more constrained singlet mixing case 
is described by the same Lagrangian with k = sin^ 9 < 0.15. In this simplified model framework we do not introduce 
a direct coupling between the SM Higgs and XX this interaction can be easily captured with giving (j) the same 
mass as the Higgs. 


III. COMMENTS ON THE RELIC DENSITY AND DIRECT DETECTION CONSTRAINTS 


Simplified models for dark matter are introduced to capture the main aspects of dark matter collider phenomenology, 
without being complete models. It is therefore customary not to impose constraints from relic density or direct 
detection stringently. Still, the model introduced above in Eq. (12) is a valid model that could have cosmologically 
viable dark matter. Therefore, to give an indication of constraints for models of this type, we calculate the relic 
density and direct detection constraints assuming that 


(i) the dark sector fermions y x which enter the simplified model definition (121 is the cosmologically stable dark 


matter and not merely one of the dark sector degrees of freedom which are long-lived on a collider scale; 


(ii) the dark matter particles annihilate predominantly via the mediator interaction specified in (121, and there 
are no other DM annihilation channels beyond the simplified model (12) or that they are highly suppressed. 


We stress that if either of these additional assumptions are not satisfied, the relic density- and direct detection-related 
constraints discussed in this section will not apply. These are strong assumptions, that can easily be evaded in many 
well motivated DM models. 

We will now require that the dark matter does not overdose the Universe, and that the direct detection cross section 
is sufficiently small to not having been observed so far. 

We calculate the relic density and direct detection limits using the MadDM [5^ with the simplified model 
( [T^ . The computed relic density is compared to the observed relic density from the Plank Satellite [57j of = 
0.1199 ± 0.0027, and the direct detection cross sections are compared to the limits from the LUX experiment [2H]. 
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Figure shows the contours of the computed relic density and the direct detection exclusions on the mediator mass 
- dark matter mass plane and for various values of 5 dm and k. For the direct detection constraint we have assumed 
that the DM density interacting with the detector is given by the canonical value, even if the DM in our model only 
is a sub-component of the total DM density in this region of parameter space. Therefore, the direct detection limits 
on our model are weaker than what is shown in the figure in the region of parameter space where the calculated DM 
density is smaller than the observed value. 

For the collider phenomenology at the LHC and at future colliders, we will be interested in heavy mediators with 
the dark matter mass and the dark matter coupling largely unconstrained as long as the scalar mediators have a 
large branching ratio to dark matter. From Figure we can see that all these models easily avoid direct detection 
constraints, and as long as the dark matter mass is quite heavy we can have mediator masses up to 2500 GeV, without 
over-closing the Universe. For 5 dm = 4 and a heavy mediator Mmed — 2.5 TeV we need > 400 GeV to have viable 
dark matter (another way to put it is that only the DM which is more than 20 times lighter than the mediator is 
constrained here). For smaller couplings the minimal DM mass increases accordingly (as can be seen from the second 
and third plots in Fig. to not overdose the universe for the heaviest mediators, but this is not a problem for the 
models we will consider in the rest of the paper. We will therefore now turn to collider phenomenology where we will 
study models that, if we interpret them as complete models, can provide a viable dark matter candidate. 

We conclude that the relic density and direct detection considerations can provide useful constraints on our simplified 
model under certain assumptions, in addition to the collider searches. This provides an important complementarity 
to the collider phenomenology we will now study. If the LHC or future colliders can resolve and probe the mediator 
mass-scale and a signal with missing energy is discovered, one of the main open questions will be if the signal results 
in the production of cosmological dark matter and what is its particle identity. 


IV. COLLIDER LIMITS ON SCALAR MEDIATORS WITH 2 JETS AND MISSING TRANSVERSE 

ENERGY AT THE LHC 


For deriving the collider limits on the models with scalar mediators to dark matter sectors and for the ability 
to distinguish between models with different mass scales, we will use the search strategy based on final states with 
missing transverse energy plus two jets. There are four main kinematic quantities associated with the plus 2 
jets signatures: the missing transverse momentum the jets invariant mass Mjj, the azimuthal angle between the 
tagging jets and the jets pseudo-rapidity difference Ary = ry^i — 5 ^ 2 • In terms of these we impose the VBF cuts 

[21 ED], 

> 100 Gev, > 1200 GeV, A()ijj < 1, Ary > 4.5, VT,j > 40 GeV, (13) 


to separate the signal and background. Here is the transverse momentum of each jet defined by using the anti-kt 
jet algorithm with R = 0.4. We reconstruct jets using Fastjet m- After imposing these cuts, the main production 
channel of the scalar mediator is largely reduced to weak vector boson fusion (WBF), leaving only a small contribution 
from the gluon fusion (GGF) channel (c/. Table 1). In spite of the relative smallness of the GGF process after cuts 
(13), one should not be tempted to approximate them by the Higgs-gluon effective vertex. The inclusion of finite top 
quark mass effects in the top-loop in the GGF production is known to be important (in the context of DM searches at 
the LHG this was emphasised in [32)1. especially for heavier scalar mediators where the heavy top mass approximation 
breaks down. We therefore simulate both the WBF and GGF contributions to the signal with VBFNLO [551 - 155] which 
includes the full top-loop dependence to GGF. 

The background is simulated at leading order using MadGraph[56]. Both signal and background are then showered 
with Herwig-|--|- m- The main backgrounds are Z -|- 2 jets with the Z decaying to neutrinos and W^ -I- 2 jets where 
the W decays to a neutrino and a missing lepton. We count the lepton as missed if it has |7yi| > 2.5 or px < 10 GeV. 
We have also checked that the it background is negligible after the cuts. The projected LHC exclusion limits for these 
final states have been studied previously in |29l EQ] EH |39] in the context of an invisible branching ratio for the SM 
Higgs. 


A. Width effect on differential distributions 

In Figure we can see the effect of varying the width of the mediator on the differential distributions of Mjj 
and A(j)jj for a mediator with Mmed = 800 GeV. A smaller width leads to a slightly broader Mjj tail and flatter 
A(j)jj distribution. For reasonably small total widths this effect is not very large. We will therefore use the narrow 
width approximation, where we produce the mediator on-shell with subsequent decay to xx with a branching ratio 
determined by the coupling constants and dark matter mass, when we simulate the signal. 
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FIG. 2: Dark matter relic density and direct detction constraints for our simplified model for dark matter for various values 
of guM and k. The lines give relic density contours and the grey region shows the area exclude by direct detection constraints. 




FIG. 3: Kinematic distributions for different values of the mediator width at yfs = 13 TeV when Mmed = 800 GeV. 
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B. The LHC exclusion limits reach 

Our first goal is to establish the projected LHC exclusions for models with scalar mediators based on the 2 jets and 
final states. We aim to evaluate the upper limit on the mediator mass for the model to be within the LHC reach. 

The left panel in Table 1 shows the cross-sections for the signal at the LHC at the 13 TeV center-of-mass energy, 
assuming a 100% branching ratio of the scalar mediators to dark matter and k = 1. The cross-sections for SM 
backgrounds are shown in the table on the right. Using these one can calculate the simple projected exclusion limits 
for these models from a standard cut-and-count procedure. 


A7med 

VBF 

GGF 

Total 

125 GeV 

89 

17 

107 

250 GeV 

61 

13 

74 

500 GeV 

26 

10 

36 

750 GeV 

12 

3 

15 

1000 GeV 

6 

0.7 

6.7 

1500 GeV 

2 

0.08 

2.01 


Background 

Gross-section(fb) 

Zjj 

128 

W+jj 

116 

w-jj 

40 


Table 1. Cross-sections (fb) at partonic level after VBF cuts in (13) at 13 TeV. 


For our analysis we will use the differential cross-sections to perform a binned log-likelihood analysis [40] to com¬ 
pute confidence levels (CLs) for experimental exclusions [41]. In the four plots of Fig. we show the normalised 
differential distributions for signal and background as functions of the four kinematic variables Mjj, Aij, A(pjj. 
These kinematic distributions are plotted for different values of the mediator mass ranging from Mj^ed = 125 GeV to 
1500 GeVb 

The differences in their shapes for models with different values of Mmed can be used for differentiating between them. 
The binned log-likelihood technique for computing confidence levels is based on regarding each bin in a histogram 
for the measured variable as an independent search channel to be combined with all others. Systematic uncertainty 
is taken into account by running many pseudo Monte Carlo experiments where the normalisation of the background 
histogram is varied randomly. The significance is then given by the fraction of these experiments that has a smaller 
likelihood ratio than that for the expected background distribution. We use the Mjj distribution(the Arj distribution 
gives similar results) with ten bins to both determine the signal exclusion limits and later to distinguish between 
different signal models. 

In Figures [5]|^ we show the LHC reach for excluding scalar mediator models for different values of mediator masses. 
Figure [^applies to generic models with k = 1 and assumes a 5% and a 10% level of systematic uncertainty. Figure]^ 
shows the LHC exclusion contours in the context of the mediator-Higgs mixing models, here we set k = 0.15 and 
assumes a 1% and a 5% systematic uncertainty. Plots in Fig. show the LHC exclusion limits without fixing the k 
parameter to a specific value. Here one allows k to float such that for each model the computed cross-section is set 
equal to a cross-section that corresponds to a 30% invisible branching ratio for the 125 GeV Higgs. 

The conclusions we draw is that for generic scalar mediator models with k ~ 1, with the 13 TeV LHC we can probe 
models with mediator masses up to Mmed ~ 750 GeV (assuming a 5% level of systematic uncertainty). For the models 
with small k, in particular the models associated with the Higgs-singlet-mediator mixing where k = sin^ 9 < 0.15, we 
can probe up to M^^d ~ 500 GeV (with an optimistic 1% systematic uncertainty). The decrease in cross-section at 
small values of k not surprisingly makes it very hard to reach to the higher mediator masses in the Higgs portal-type 
mixing model realisations at the LHC. While we only focus on CP-even scalar mediators, the same distributions can 
be used to access the scale of CP-odd scalar mediators. A(j)jj can be a powerful observable to discriminate CP-even 
from CP-odd mediators (UHlll- 


C. Distinguishing between models with different mediator masses 

For the models which are within the LHC reach, i.e. with Mmed below the upper bounds set be the exclusion 
contours in Figs. [5]|^ the next step is to be able to distinguish between different models. 


^Compared to differential distributions at e+e colliders mils], at the LHC differences between the models are less pronounced and more 
difficult to exploit. 
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FIG. 4: Kinematic distributions for different values of the mediator mass for the signal, and for the background at the LHC. 
Mjj distributions are shown on the top left panel, is on top right right, A 77 and (f>jj distributions are on the bottom left and 
right panels respectively. 




FIG. 5: We characterise the LHC reach for models with different values of M^ed by computing confidence levels for excluding 
signals from the SM backgrounds. We consider models with k = 1 and on the left panel use a systematic uncertainty of 5%, 
the panel on the right corresponds to 10 % systematic uncertainty. 
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FIG. 6: The LHC reach for different Mmed models with k = 0.15 in terms of confidence levels to exclude signal from SM 
background. On the left panel we use a systematic uncertainty of 1%, and the panel on the right corresponds to 5% systematic 
uncertainty. 




FIG. 7: The LHC reach for different Mmed models. We normalise the cross-sections for all models to the SM Higgs cross-section 
with Brinv = 30% and a systematic uncertainty of 5% (left panel) and 10% (right panel). 


This is achieved by comparing the shapes of the kinematic distributions plotted in Fig. for different mediator 
masses. As we increase Mmed the visible jets will recoil against a heavier object and we can use the kinematic variables 
to distinguish the models. We will use the following procedure. Before we even start comparing different models, we 
would need an excess of signal events over the SM background in the data after the VBF cuts. The cross-section of 
this signal can be used to infer an upper limit for the mass of the mediator as a function of k. The question then 
becomes if we can distinguish between the different models that can achieve the measured cross-section. We will again 
use a binned log-likelihood method and we will be comparing models pairwise. For each cross-section we select the 
two extreme models: first one with the maximal mass, and the second (reference model) - with the 125 GeV mediator. 

In Figure we can see how well one can differentiate the models at the LHC with k = 1. Specifically, all the models 
with Mmed = 250, 500 and 750 GeV can be distinguished from the 125 GeV mediator. Within our approach this 
conclusion is valid even with a relatively high systematic error of 10%. We also note that for k = 0.15 it is no longer 
possible to differentiate any of the models since the cross-section becomes too small. 

So far in Fig. we have characterised the simplified model signals by fixing the scaling parameter k to either 1 or 
0.15. Alternatively we can set the signal cross-section to a fixed value corresponding to a 30% invisible branching ratio 
of the SM Higgs. This is shown in Fig.|^which leads one to conclude that the models with heavier and heavier mediator 
masses are easier and easier to distinguish from the reference model. When comparing models with different mediator 
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FIG. 8: Differentiating the models at k = 1 at the LHC. For each valne of M^ed between 250 and 1500 GeV the models are 
compared to the reference model with the 125 GeV mediator. We assume a systematic uncertainty of 5% (left panel) and 10% 
(right panel). 




FIG. 9: Differentiating the models with the floating k parameter defined as in the caption of Fig. [7] Models are compared 
pairwise to the 125 GeV reference model. 


masses, there are in general two competing effects: the increased difference in the shape of differential distributions, 
and the decrease in the cross-section with the increase of the mediator mass. By fixing the cross-sections in Fig. the 
differences between the models are only due to the shapes of differential distributions, while in Fig.j^both effects are 
important. This explains why for example the model with Mmed = 500 GeV is easier to distinguish than the 250 and 
750 GeV models. 

We also compare models where the reference model is not the 125 GeV Higgs. The results for using 250 GeV and 
500 GeV as reference models are shown in Fig. 10 In the same way as for the 125 GeV Higgs, the cross-section for 


the reference model is set equal to that of the model we compare it with. We see that the 500 and 750 GeV models 
can be distinguished from the 250 GeV model at the LHC. At the same time, the 750 GeV model (and above) cannot 
be distinguished from the 500 GeV reference point. 
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FIG. 10: Differentiating the k = 1 models at the LHC. In the left panel we compare to the reference model with M^ed = 250 
GeV, and on the right the reference model is 500 GeV. We assume a systematic uncertainty of 5% . 


V. SCALAR MEDIATOR MODELS AT 100 TEV 


We use a similar approach to investigate the models reach and the ability to distinguish between different models 
at a future 100 TeV circular proton-proton collider. The signal and background are simulated in the same way as for 
the LHC analysis and we use the same binned-log likelihood analysis for exclusion and differentiation of the various 
models. The main difference is that we use the cuts 


> 100 Gev, Mjj > 1200 GeV, < 0-5 j A77 > 5.5, ptj > 110 GeV, 


(14) 


instead of the normal VBF cuts in (131 as we need to reduce the background more. We also allow for larger jets by 
using the anti-kt jet algorithm with R = 0.8. 

The left panel in Table 2 shows the cross-sections for the signal at a future 100 TeV collider, assuming a 100% 
branching ratio of the scalar mediators to dark matter. 


fklmed 

VBF 

GGF 

Total 



125 GeV 

120 

100 

220 



500 GeV 

76 

no 

185 

Background 

Gross-section(fb) 

750 GeV 

53 

46 

99 

Zjj 

239 

1000 GeV 

40 

20 

60 

W+jj 

298 

1500 GeV 

22 

5 

27 

w-jj 

157 

2000 GeV 

14 

2 

16 



2500 GeV 

8.7 

0.6 

9.4 




Table 2. Cross-sections (fb) at partonic level after the cuts (14) at 100 TeV. 


We have set k = 1 and assumed 100% branching ratio of the mediator to XX DM. The cross-sections for SM back¬ 
grounds are shown in the table on the right. 


Figure [TT| plots the kinematic distributions for Mjj, Ary, for models with different values of Mmed and the 
SM background. We first investigate the exclusion limits at 100 TeV for invisible decays of the 125 GeV Higgs. In 
Figure [T^ we plot exclusion limits on the branching rations to invisibles and conclude that for small systematic errors 
one can exclude Brinv ^ 10%. 

We can now consider simplified models with different values of the mediator mass. In Figure we show the 
expected reach for a 100 TeV collider for these models with rt = 1 for 1% and 5% level of systematic uncertainty. 
Figure [14] gives the expected exclusion limits for models with k = 0.15. From these Figures we conclude that in the 
case of small systematic uncertainties the 100 TeV collider would provide a very significant increase in the exclusion 
reach for new physics models with mediator masses up to 2.5 TeV. For a larger systematic uncertainty at the level of 
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FIG. 11: Kinematic distributions for different values of the mediator mass for the signal, and for the background at a 100 
TeV collider. Mjj distributions are shown on the top left panel, is on top right right, Arj and (pjj distributions are on the 
bottom left and right panels respectively. 



FIG. 12: 100 TeV reach for excluding invisible decays of the 125 GeV Higgs boson. On the left panel we use a systematic 

uncertainty of 1%, and the panel on the right corresponds to 5% systematic uncertainty 
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FIG. 13: 100 TeV reach for different Mmed models with n = 1. On the left panel we use a systematic uncertainty of 1%, and 

the panel on the right corresponds to 5% systematic uncertainty. 




to 

O 


10-2 



Systematic : 0.05, k = 0.15 


- filmed = 500 GeV 

- M,ned = 750 GeV 

- Aimed = 1000 GeV 

- = 1500 GeV 

- M,ned = 2000 GeV 

- A/,ned = 2500 GeV 


lO-^L 


20 40 60 80 


Luminosity [fb 


100 


FIG. 14: 100 TeV reach for different Aimed models with k = 0.15. On the left panel we use a systematic uncertainty of 1%, 

and the panel on the right corresponds to 5% systematic uncertainty. 


5% the reach in Aimed is 1 TeV. For the more restricted case of k = 0.15 Higgs-mixing models, the FCC could exclude 
models with up to 750 GeV. 


Finally, as we have done at the LHC energies before, we can use the kinematic distributions at 100 TeV to differen¬ 
tiate pairwise between models with different mediator masses. We first compare each model with a reference model 



reference models corresponding to mediator masses of 500, 750, 1000 and 1500 GeV. We conclude that for all the 
reference models with up to 1 TeV mediator masses we can successfully distinguish all higher mass reachable models. 
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FIG. 15: Differentiating the models at ac = 1 at 100 TeV. For each value of Mmed between 500 and 2500 GeV the models are 
compared to the reference model with the 125 GeV mediator. We assume a systematic uncertainty of 1% (left panel) and 5% 
(right panel). 




FIG. 16: Differentiating the models at k = 0.15 at 100 TeV. For each valne of Mmed between 500 and 2500 GeV the models 
are compared to the reference model with the 125 GeV mediator. We assume a systematic uncertainty of 1% (left panel) and 
5% (right panel). 


VI. SUMMARY AND CONCLUSIONS 

We have studied collider limits for simplified models of dark matter. In the case we consider, dark matter is coupled 
to the Standard Model by a scalar field mediator which interacts with both, the SM fermions and the vector bosons. 
Our main goal was two-fold, first to determine the projected reach of the Run 2 LHC and the future circular hadron 
collider for excluding such models. Second, was to use the kinematic distributions to be able to extract information 
on the values of the mediator masses. To achieve this we focussed on the two jets plus missing transverse energy hnal 
states. In our case both production mechanisms are important: the gluon fusion process which effectively probes the 
coupling of mediators to quarks, and the weak vector boson fusion which is due to the interactions with vector bosons. 

We found that at the 13 TeV LHC one can probe dark matter simplihed models with mediator masses up to 750 
GeV. At a 100 TeV collider the reach is increased to 2.5 TeV mediators. These measurements are fairly dominated 
by systematic uncertainties which we have kept relatively small and have varied between the 1% and 10% level. 

In order to ‘measure’ the mediator mass - which is the most relevant new physics scale of these models at colliders 
- we have compared the models with different values of Mmed pairwise. In both cases, at 13 and at 100 TeV, we 
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FIG. 17: Differentiating the models at k = 1 at a future 100 TeV Collider. For each value of Mmed between the reference 
model mass and 2500 GeV the models are compared to the reference model with the 500 GeV mediator in the left panel and 
750 GeV in the right panel. We assume a systematic uncertainty of 5% . 




FIG. 18: Differentiating the models at k = 1 at a future 100 TeV Collider. For each value of Mmed between the reference 
model mass and 2500 GeV the models are compared to the reference model with the 1000 GeV mediator in the left panel and 
1500 GeV in the right panel. We assume a systematic uncertainty of 5% . 


found that we can differentiate between essentially all the discoverable benchmark models we have considered. In 
particular, we can distinguish a 125 GeV reference model from models with mediator masses of 250, 500 and 750 GeV 
at the LHC. For the 100 TeV case the list of distinguishable models extends to 2.5 TeV. 

For the more constrained case of the Higgs portal models, where scalar mediators mix with the SM Higgs and the 
mixing angle is small, sin^ 6 = 0.15, the collider reach is reduced due to the overall decrease in the signal cross-section. 
We found that the LHC reach is 500 GeV for an optimistic 1% systematic uncertainty, and this is increased to 750 
GeV at a 100 TeV collider. 
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